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Abstract. Thermal conditions in the LTP, the LISA Technology Package, 
are required to be very stable, and in such environment precision temperature 
measurements are also required for various diagnostics objectives. A sensitive 
temperature gauging system for the LTP is being developed at lEEC, which 
includes a set of thermistors and associated electronics. In this paper we discuss 
the derived requirements applying to the temperature sensing system, and address 
the problem of how to create in the laboratory a thermally quiet environment, 
suitable to perform meaningful on-ground tests of the system. The concept is a 
two layer spherical body, with a central aluminium core for sensor implantation 
surrounded by a layer of polyurcthane. We construct the insulator transfer 
function, which relates the temperature at the core with the laboratory ambient 
temperature, and evaluate the losses caused by heat leakage through connecting 
wires. The results of the analysis indicate that, in spite of the very demanding 
stability conditions, a sphere of outer diameter of the order one metro is sufficient. 
We provide experimental evidence confirming the model predictions. 
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1. Introduction 



LISA Pathfinder (LPF) is an ESA mission, with NASA contributions, whose main 
objective is to put to test critical parts of LISA (Laser Interferometer Space Antenna) , 
the first space borne gravitational wave (GW) observatory The science module 
on board LPF is the LISA Technology Package (LTP) [2] , which basically consists in 
two test masses in nominally perfect geodesic motion (free fall), and a laser metrology 
system; this one detects residual deviations of the test masses' actual motion from the 
ideal free fall, to a given level of accuracy 

In order to ensure that the test masses are not deviated from their geodesic 
trajectories by external (non-gravitational) agents, a so called Gravitational Reference 
System (GRS) is used This consists in position sensors for the masses which send 
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signals to a set of micro-thrustcrs; the latter take care of correcting as necessary the 
spacecraft trajectory, so that at least one of the test masses remains centred relative to 
the spacecraft at all times. The combination of the GRS plus the actuators is known 
as drag-free subsystcm||. 

The drag-free is of course a central component of LISA, and needs to be operated 
at extremely demanding levels of accuracy. The laser metrology system should then 
be sufficiently precise to measure relative test mass deviations. The overall level of 
noise acceptable for LISA is defined in terms of rms acceleration spectral density, and 
has been set to PP 



<lisaM< 3x10-1^ 



1 



3 mHz / 



ms-^/VlE (1) 



1 /2 

in the frequency range 0.1 mHz < uj/2tt < 100 mHz. This is equivalent to 5^ ~4 x 
10-21 Hz-i/^ with the same frequency dependence. 

LPF is conceived, as mentioned above, as an in-flight test of key technologies for 
LISA. Top level performance requirements for LPF have however been relaxed by an 
order of magnitude relative to LISA, both in noise amplitude and in frequency band, 
to still challenging goals [5]: 



^yLPFM< 3x10-14 



1 + 



Lu /2Tr Y 
3 mHz J 



ms-^/Vu^ (2) 



in the frequency range 1 mHz < lu/2tt < 30 mHz. We shall be referring to this 
frequency band as the LTP Measuring Bandwidth (MBW) in the sequel. 

Equation (01 gives the global noise budget. This is made up of contributions from 
various perturbative agents, both of instrumental and environmental origin. One of 
these is temperature fluctuations, for which a stability requirement has been set to 

S}/^ {uj) <10-*K/Vb^ , ImHz < w/27r < 30mHz (3) 

in order to comply with |(5Jl under suitable noise apportioning criteria — see section |21 
Because temperature stability is important, a decision has been taken to place 
high precision thermometers in several strategic spots across the LTP — as part of what 
is called Diagnostics Subsystem Such high precision temperature measurements 
will be useful to identify the fraction of the total system noise which is due to thermal 
fluctuations only, and this will in turn provide important debugging information to 
assess the performance of the LTP [7] . 

The best temperature sensors for our purposes are electric devices whose ohmic 
resistance varies with temperature. We have chosen to use thermistors, which are 
semiconductor resistors whose resistance decreases as their temperature increases 
because they have a rather steep sensitivity curve, and this suits our needs. Such 
sensors need additional electronic circuitry to bias them and acquire data. The entire 
chain of sensors plus electronics must of course be tested in ground before boarding, 
and this requires suitably stable environmental conditions in the first place. 

This paper addresses the problem of which are these conditions, and how to 
implement them for a reliable laboratory test of the LTP temperature sensors and 

II The term drag-free dates back to the early days of space navigation, when it was used to name a 
trajectory correction system designed to compensate for the effect of atmospheric drag on sateUites 
in low altitude orbits. 

^ The Diagnostics Subsystem of the LTP also includes magnetometric measurements and a charged 
particle flux detector. 
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electronics. It is organised as follows: in section|21we review and quantify the identified 
sources of thermal noise which result in the budget given in equation (j2J). In section 13 
we discuss how the stability requirement results in a requirement on the temperature 
system performance, and extend the argument to make precise the environmental 
conditions which must be met in the laboratory test. In section 0] we present the 
physical hypotheses, lay down the mathematical model of the proposed insulation 
scheme, and find an analytic solution to the equations. We then address in section 
the numerical implications of the above for realistic situations, including a discussion 
of heat leakage through connecting wires. Section |H1 deals with the experimental 
verification of the model predictions, and section [T] summarises our conclusions and 
future prospects. Some supplementary mathematical detail is provided in an appendix. 

2. Thermal disturbances in the LTP 

Temperature fluctuations inside the LTP result in noisy readout at the interferometer 
output port — the phasemeter. The reason for this is that there are components in 
the LTP whose behaviour is sensitive to temperature changes, as we shall describe in 
this section. 

As a rule of thumb, the total contribution of thermal noise to the total acceleration 
noise, equation ||2Jl, should not exceed 10%. We thus require that 



for frequencies within the LTP MBW. This assumption is in fact somewhat 
conservative, as the Project Engineers have estimated that more than twice this value 
is actually compliant with the overall LTP noise budget [H]. We shall however adopt 
equation Q as reference to ensure we are playing on the safe side. 

For the sake of clarity, we consider separately the influence of temperature on the 
GRS and on the Optical Metrology System (OMS). 

2.L Noise effects inside the GRS 

Temperature differences between the walls of the electrode housing cause differential 
pressures on opposite faces of the test masses, which in turn result in net forces on 
them, hence in noise at the phasemeter. 

Three different mechanisms have been identified whereby temperature fluctua- 
tions distort the LTP readout: radiation pressure, radiometer effect and outgassing. 
Let us briefly describe each of them, and quantitatively estimate their respective con- 
tributions to thermal noise. 

2. LI. Radiation pressure A body at any (absolute) temperature T emits thermal 
radiation. This exerts pressure on any surfaces the radiation hits. According to 
standard electromagnetic theory, such pressure is given by 



where a = 5.67x 10 * Wm ^ is the Stefan-Boltzmann constant, and c is the speed 
of light. Thus, if there are temperature fluctuations around the test mass, a noisy net 




Pq . m . 




(5) 
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Electrode housing 



Figure 1. Schematics of tlie effect of different pressures on opoosite faces of a 
test mass. The physical origin of pressure differences is explained in the text. 



force will appear on it — see Figure for a graphical display. The effect can be easily 
quantified making use of equation jS)): 

Ape..„. = ^T3Ar (6) 
3c 

where Ap and AT make reference to differences of pressure and temperature between 
the sides of the test mass. Associated acceleration noise is hence obtained multiplying 
the above by the test mass surface area, f^jyj and dividing by its mass totm^ 



16 ^xm'- 
3 ttt-tmc 



(7) 



2.1.2. Radiometer effect This is an effect which happens in rarefied gas atmospheres, 
its name historically coming from its association with the theory of Crookes's lightmill 
radiometer. In low pressure atmospheres, where the gas particles have a mean free 
path well in excess of the dimensions of the containing vessel, equilibrium conditions 
do not happen when pressure is uniform, but rather when the ratios of pressure to 
square root of temperature equal one another. Or 
PI P2 



The pressure gradient is readily obtained from this expression, and thence the 
associated test mass acceleration: 

2 mxM T 

2.1.3. Outgassing Outgassing is one of the causes of the presence of gas within the 
walls of the GRS. In the present contetx, outgassing problems actually derive from 
temporal fluctuations in its rate, which once more result in pressure fluctuations, 
thence in noise. Outgassing rates are very strongly dependent on materials and 
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geometry, and the theoretical analysis of the problem is by no means a simple one. 
The issue has also attracted attention of people from other space missions and 
further experimental work appears to be necessary to reliably assess the impact of 
this phenomenon. Partial evidence has however been gathered that outgassing might 
be in practice a small effect in the LTP ^^l- We shall therefore omit any further 
consideration of this phenomenon here. 

2.1.4- Total thermal noise in the GRS If we make the assumption that radiometer 
and radiation pressure fluctuations are uncorrelated, and neglect outgassing, then 
equations {Tj) and ^ are added quadratically, and hence the spectral densities of 
acceleration and temperature in the GRS arc related by 

^•s.™- ™s(") - /(iiH^(i£^ (10, 

Nominal conditions in the LTP are the following: 

= 4.6 X 10-2 m 
"T-TM = 1-96 kg 
T = 293 K 

p= 10-5 Pa 

which give 

St%bA^) = 3.3 X 10^° 5;;/L_, GRs(^) KHz-i/2 (11) 

This expression gives 0.9 x 10~^K/v^Hz in the worst case that all the thermal 
acceleration budget, equation Q), is allocated to thermal fluctuations in the GRS. 

2.2. Noise effects inside the QMS 

The Optical Metrology System of the LTP is based on a non-polarising, heterodyne 
Mach-Zender interferometer ^Sl- The light source is an infrared Nd:YAG laser 
of 1.064 ^m wavelength and approximately 0.25 W of power. Optical components 
are mounted on an optical bench of highly stable optical properties. Temperature 
fluctuations basically affect the optical system through two distinct effects |14| : 

• The index of refraction of optical components depends on temperature. 

• Temperature changes cause dilatations (and contractions) of optical elements, 
which in turn cause light's optical path to change accordingly. 

While it is not difficult to characterise how individual components are influenced 
by the above effects, to assess the behaviour of the fully integrated optical metrology 
is a more complicated task. Significant progress has been made since the early 
design proposals — see reference |2] and following articles in that issue of CQG — , 
and improved materials and designs are now available. Altogether, it appears that 

4(lMS^10''K/Vlfe (12) 

is a sensible requirement which should comfortably guarantee the performance of the 
optical bench against temperature fluctuations in flight — see chapter 12, and )15| . 
Like before, the noise level l|12(l is estimated to account for about 10% of the total 
LTP acceleration noise, equation 
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2.3. Total thermal noise budget for the LTP 

If we make the hypothesis that noise in the OMS is uncorrelated with noise in the 
GRS then the respective spectral densities add up, i.e., 

5't(w) = St, grs{^) + St, OMsi^^) (13) 

Inserting here equations Hll() and (|12|l . we find that the LTP thermal stability 
must be 

S^^^iuj) ~ 1.3 X 10"'* K/Vlfe (14) 

where the usual frequency dependence of equation has been dropped, as it 
is in practice the case that temperature fluctuations actually drop towards higher 
frequencies — due to efficient thermal insulation of the LTP. 

The number just obtained is very close to that of equation (j2J), and derived under 
worst hypotheses in each case. We are thus reassured that equation (jj)) is a sensible 
requirement for the temperature fluctuations which can be tolerated in the LTP. Let 
us however stress that we still can count with some margin, as acceleration thermal 
noise allocation has been taken as a conservative 10 % of the total acceleration noise 
— see paragraph following equation 10} . 

3. Temperature measurement and on-ground test requirements 

Temperature stability requirements for the LTP are now established. We rewrite them 
again: 

s!/'^{uj) <lO-'^K/VTh. , ImHz < w/27r < 30mHz (15) 

Because this is a requirement, satellite and payload will be built such that thermal 
conditions in the LTP meet ()15|l . Even so, to monitor the magnitude of anyway 
compliant temperature fluctuations is a powerful diagnostic tool, yielding valuable 
information on the system thermal behaviour, which will be useful for LISA. In 
addition, we wish to be able to diagnose whether the conditions H15|l actually prevail 
at any given time during the mission. For both objectives, temperature gauges are 
obviously needed. 

We now pose the question: which is the level of noise we can accept in the 
temperature measurement system — which includes sensors, wires and electronics — if 
we are to diagnose temperature variations below the level H15() ? Clearly, the answer to 
that question depends on how accurately those fluctuations arc to be measured. The 
stability requirement equation (|15|l is already rather demanding of itself, so we do not 
expect thermal fluctuations to be much smaller. With this in mind, to request the 
measuring system to be about one order of magnitude less noisy than the maximum 
noise level to measure seems a sensible option: 

4(Lasuremc„t(^) < 10^' ^/VWz , 1 mHz < cj/27T < 30 mHz (16) 
This has in fact become a mission top level requirement — see [5], section 6.2. 
There are two groups of reasons which support it: 

(i) Equation (|15|l dcflnes the maximum acceptable level of temperature noise in the 
LTP. If this is satisfled, which of course must, then actual fluctuations will be less 
than that. Requirement H16|) then sets a 10 % minimum discrimination capability 
for the measuring device, a standard approach which is certainly compatible with 
better performance. 
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(ii) LISA is more demanding than LPF as regards thermal stabihty. Actually, LISA 
requires an order of magnitude less thermal noise than LTP ^1]. If we require (|16|l 
for LTP then we are in a position where analysis of thermal sources of noise of 
relevance for LISA can be identified and tagged for improvement, at least in the 
overlapping frequency band of both missions. This prospect is in line with the 
very concept of LPF as a precursor mission. 

3.L Absolute versus differential temperature measurements inflight 

The above arguments, and specially the first, can perhaps be critisised in terms of: 
why not perform differential measurements? This might relax the very demanding 
requirement in equation H16|l . in the sense that it would only apply to differential 
rather than absolute temperature measurements in fiight. 

While it is true that certain thermal disturbances depend on temperature 
gradients across the test masses — like the radiometer effect and radiation pressure 
gradients — there are others which do not — mostly those related to the optics. One 
could accordingly split up the temperature gauges into two classes, but this does not 
seem a particularly sensible choice, since the best device would obviously be the one 
to use in all cases, anyway. 

A space mission like LPF does not generally allow to fix hardware design 
inefficiencies once it has been launched. The choice of making applicable the 
requirement stated in equation H16() to all temperature measurements, whether 
differential or absolute, seems thus not advisable to relax: some margin is necessary 
to cope with unforseen sources of error. 

3. 2. Environmental conditions for on-ground tests 

Before launch, the temperature diagnostics hardware must be tested in ground. In 
order to do a meaningful test, a sufficiently stable thermal environment must be 
granted for the sensors in the first place. Here, sufficiently stable means that any 
observed fluctuations in the readout data should be attributable solely to sensor 
and/or electronics noise, rather than to a combination of the latter with ambient 
temperature fluctuations. This means ambient temperature fluctuations in the 
thermometers should be distinctly smaller than the limit set by the the requirement, 
equation H16|l. Again, one order of magnitude below that target sensitivity seems a 
good prescription. So we require 

^'r^testbcdl'^) ^ lO^^K/Vlfe , ImHz < w/27r < 30mHz (17) 

It turns out that IO^^K/VHz is a very demanding temperature stability, orders 
of magnitude beyond the capabilities of thermally regulated rooms. Once more, this 
figure could be widely relaxed if tests were done differentially, i.e., between pairs of 
sensors close enough to one another. While such differential measurements are indeed 
envisaged, we also want to make sure that the sensing system is sensitive enough that 
absolute temperature measurements are compliant with the most exigent requirements. 
We thus adopt equation 1)17(1 as the design goal baseline. 

The concept idea of the insulator is displayed in figure |21 an interior metal core 
of good thermal conductivity is surrounded by a thick layer of a poorly conductive 
material. The inner block ensures thermal stability of the sensors attached to it, 
while the surrounding substrate efficiently shields them from external temperature 




Figure 2. Thermal insulator design concept. Left: 3D diagramme, including 
sensor placement principle. Right: cut view, with notation convention dictionary. 



fluctuations in the laboratory ambient. We propose a spherical shape for the sake of 
simplicity of the mathematical analysis, even though this will be eventually changed 
to cubic in the actual experimental device for practical reasons. 

4. Mathematical model of the insulator 

The basic assumption of the mathematical analysis we shall present is that heat flows 
from the interior of the insulator to the air outside, and from the latter to the interior 
of the insulator, only by thermal conduction. This is a very realistic hypothesis in 
the context of the experiment, as radiation mechanisms are certainly negligible and 
convection should not play any significant role, either, since the entire body is solid, 
and temperature fluctuations will be small anyway. 

Let then T(x, i) be the temperature at time t of a point positioned at vector x 
relative to the centre of the sphere. Under the conduction-only hypothesis, T(x, t) 
satisfies Fourier's partial differential equation 

pcp^r(x,t) - V- [Kvr(x,i)] (18) 

where p, Cp and k are the density, specific heat and thermal conductivity, respectively, 
of the substrate. We shall assume these are uniform values within each of the two 
materials making up the insulating body, with abrupt changes in the interface. We 
can thus represent them as discontinuous functions of the radial coordinate, as follows: 

p c kU)-[ /'I'^Pi''^! if 0<r<ai 
p,Cp,^W_ I p,^cp2,n2 if a,<r<a2 

with r= |x|. Initial and boundary conditions are the following: 

r(x,i = 0) = 0, T{r = a2,t)^Toi9,^;t) (20) 

where 9 and (p are spherical angles which define positions on the sphere's surface. The 
boundary temperature can be expediently expressed as a multipole expansion: 

To(0,¥>;t) =^ &z,„(i)ri,„(0,<^) (21) 

Ini 
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where Yi^n^{9,Lp) are spherical harmonics, and bi„^{t) are boundary muhipole 
temperature components. 

In practice, the boundary temperature will be randomly fluctuating, therefore 
birnit) will be considered stochastic functions of time. We shall also reasonably assume 
them to be stationary Gaussian noise processes with known spectral densities, Sira{^)- 

As shown in the appendix, the frequency analysis of this problem leads to a 
transfer function expression of the temperature inside the body: 

T(X, = ^ Hlra{^, oj) biM (22) 
Im 

where tildes {") stand for Fourier transforms, e.g., 

/oo 
r(x, t) e-"^' dt (23) 
-oc 

etc. If we make the further assumption that different multipole temperature 
fluctuations at the boundary are uncorrelated, i.e., 

(J)*l,^,{uj) bl„i{uj)) = Slrn{u}) Sl>l Srn'ni (24) 

then the spectral density of fluctuations at any given point inside the insulating body 
is given by 

5't(x,w) = ^ \Hi,n{x,u!)f Si,n{co) (25) 

Irn 

It is ultimately the spectral density ^^(x, oj), for |x| <ai, which has to comply 
with the requirement expressed by equation 1)17(1 . Based on knowledge (by direct 
measurement) of ambient laboratory temperature fluctuations, equation ((25(1 will 
provide the guidelines, as regards materials and dimensions, for the actual design 
of a suitable insulator jig. 

Isotropic boundary conditions 

Thermal conditions in the laboratory are rather homogeneous. This means that the 
boundary temperature fluctuations will be in practice essentially independent of the 
angles and (p, i.e., 

To{e,cp;t)^Bit) (26) 

and consequently the generic expansion equation ((21(1 includes only the monopole 
term. Hence 

6oo(t) - VI^B{t) (27) 

The temperature r(x, cj) in this case will only depend on radial depth, r, and 
therefore 

f{r,Lj) = H{r,uj)B{uj) (28) 
with iJ(r, tj) = VItt _ffoo(x, w). According to equation l(A.21() of the Appendix, this is 

i/(.,.)^|^"^"^^°^"^^^' "-^-"^ (29) 

[ rio{uj) ja{j2r) + Co(w) Vo{l2r) , ai < r < 02 

This is a low-pass filter transfer function — even though the cumbersome 
frequency dependencies involved in the expressions above do not make it immediately 
obvious. The 3dB cut-off angular frequency for this filter defines a time constant t 
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Figure 3. Frequency response at the centre (r = 0) of a spherical thermal 
insulator — as given in equation 1291 — , along with the frequency response of 
a low-pass filter of order one and same frequency cut-off, a;t.ut-off = !/''"■ 



by 1/r = a;cut-off, which correspondingly is a comphcatcd function of the insulator's 
physical and geometrical properties. Figure 13 shows graphically the situation: the red 
curve plots the square modulus of H{r = 0,uj), where the cut-off l/r is shown as the 
3 dB abscissa. For comparison, the figure also shows (blue curve) a low-pass filter of 
the first order with the same frequency cut-off, i.e., |iJist ordcr(w)P = (1 + w^r^)^^. 

The high frequency roll-off of the real filter (in red) is seen to drop below the first 
order counterpart (in blue): the latter clearly has a slow slope at high frequencies, 
l^^ist ordcr(w)| ^ (wt)~^, while the former can be shown to follow a much steeper, 
exponential curve: 



As already mentioned in section^] to test the temperature sensors and electronics 
we need a very strong noise suppression factor in the LTP frequency band. Inspection 
of figure and equation l|Sn|) readily shows that high damping factors require such 
frequency band to lie in the filter's tails. The thermal insulator should therefore be 
designed in such a way that its time constant r be sufficiently large to ensure that 
the LTP MBW frequencies are high enough compared to (27rT)~^. The exponential 
roll-off in the transfer function shown by H30|l makes the filter actually feasible with 
moderate dimensions. 

5. Numerical analysis 

In this section we consider the application of the above formalism to obtain useful 
numbers for the implementation of a real insulator device which complies with the 
needs of the experiment. 

First of all, a selection of an aluminium core surrounded by a layer of polyurethane 
was made. Aluminium is a good heat conductor and is easy to work with in the 
laboratory; polyurethane is a good insulator and is also convenient to handle, as 
it can be foamed to any desired shape from canned liquid. Other alternatives are 
certainly possible, but this appears sufficiently good and we shall therefore stick to 





(30) 




this specific choice. The relevant physical properties of aluminium and polyurethane 
are as follows: 



Density Specific heat Thermal conductivity 

p (kgm-3) Cp (Jkg-iK-i) K (Wm-iK-i) 

Aluminium 2700 900 250 

Polyurethane 35 1000 0.04 



FigureElplots the amplitude damping coefficient of the insulator block, |i7(r, cj)!, 
at the lower end of the LTP frequency band (i.e., 1 mHz) at the interface position, 
r = ai. Each of the curves corresponds to a fixed value of the latter, and is represented 
as a function of the outer radius of the insulator. This choice is useful because the 
sensors are implanted for test on the surface of the aluminium core, and also because 
at higher frequencies thermal damping is stronger. The plot clearly shows that the 
assymptotic regime of equation (|30|) is quite early established. 

The choice of dimensions for the insulating body must of course ensure that the 
minimum requirement, equation H17() is met. For this, a primary consideration is the 
size of the ambient temperature fluctuations in the site where the experiment is made. 
Dedicated measurements in our laboratory showed that 

'S't! room('^) - 10"^ K/Vlfe , 1 mHz < tj/27r < 30 mHz (31) 

We therefore need to implement a device such that \H{ai,uj) \ < 10~^ throughout 
the MBW. Suitable dimensions can then be readily read off figure 0] and various 
alternatives are possible, as seen. Before making a decission, however, we need to 
make an additional estimate of the heat leakage down the electric wires which connect 
the temperature sensors with the electronics, which lies of course outside the insulator. 
We come to this next. 
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5.1. Heat leakage through connecting wires 

We use a simple model, consisting in assuming the connecting wires behave as straight 
thin rods which connect the central aluminium core with the electronics, placed in 
the external laboratory ambient. Because the polyurethane provides a very stable 
insulation, we can neglect the lateral flux, hence only a unidimensional heat flow 
needs to be considered. For this, the following equation relates the heat flux to the 
temperature difference between the two wires' edges: 

Qit) = Kwirc [T{a2,t) ~ T(ai, t)] (32) 

where Kwiro is the thermal conductivity of the wire, i?wirc its transverse radius, and 
£wirc its length inside the polyurethane layer. 

On the other hand, the heat flux results in temperature variations in the metal 
core, given by 

dT 

Q{t)^ piCpiVi—{ai,t) (33) 

where Vi —Anai/S is the volume of the metal core. Equating the above expressions 
wc find 

Atwhe^^^ [r(a2,<)-T(ai,t)] =piCpiFi^(ai,i) (34) 

For fluctuating temperatures, we can now obtain the relationship between the 
spectral density at the aluminium core and the ambient, caused by heat conduction 
along the wire: 

5't( wiro(ai' = l^wirc(w)| 5'y(^^,„bicnt (^) (35) 

where 

\H^U-)\ ^ "wi.ei?Le (36) 
LU piCpiVi £wire 

and where the approximation has been made that the temperature fluctuations at the 
inner end of the wire are much smaller than those at the outer end, due to the presence 
of the polyurethane layer. 

In practice, there will be several sensors for test inside the insulator. Under the 
hypothesis made that no lateral heat flux is relevant, the transfer function for a bundle 
of N of wires is, at most, N times that of a single wire. Thus, 

Iff 1^, ,M — ^wirc Rwhc (or7\ 
\HNwircs((^)\ - —TTT- ^ —ff (.3') 

Let us consider numerical values in this expression. We use thin copper wires 
(recu = 401 Wm^^K^^) of radius i?wirc = 0.1 mm, and assume some fiducial parameters 
for the size of the aluminium core, ai, the wire length, ^wirc, the number of connecting 
wires, N, and the frequency, iojl-K. The following obtains: 



\B^^:u^)\ = 1.1 X 10-^ ( ^\ ( V ( 

' 1^ V30/V13cm/ V25cmy V 1 mHz / 



(38) 



This result indicates that, for laboratory fluctuations in the level of equation 
leakage through wiring causes fluctuations in the temperature sensors of about 
10~^ K/VHz, equation (|35|l . which is compliant with the requirement of stability, 
equation l)17|l . The most sensitive parameter in the above expression is the size of the 
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metal core, and this determines the need to make it somewhat large. The length of 
the wires has been taken to be 25 cm, but this does not necessarily mean we need 
02 = 38 cm (assuming the radius of the aluminium core is ai ~ 13 cm), because the 
wires can be partly wound inside the polyurethane layer to further protect the system 
against leakage. In fact, wire lengthening is an easy way to improve attenuation. 

As regards frequency dependence, compliance is guaranteed in the entire MBW 
if it is at its lower end: indeed, not only |iJwiro(^)| decreases as iu~^ — see 
equation H38|l — , also ambient noise fluctuations happen to drop below 10^^K/\/Hz 
at higher frequencies. 

6. Experimental verification 

Direct experimental verification of the predicted performance of the insulator would 
require the sensing system, i.e., sensors plus electronics, to have a level of noise below 
10~^K/-\/Hz in the MBW. The requirement on the latter is however 10~'^K/\/Hz, 
and this is itself to be put to test. 

We shall nevertheless argue that the thermal environment of the sensors produced 
by the insulator is in fact better than required. Evidence of this comes from rounds 
of measurements made over periods of several days to weeks, as we now describe. 

Already processed data are displayed in figure |31 The noisy curves are estimates 
of the insulator transfer function as derived from two different experimental runs: 




10"" 10 ' 10'" 10 = lO-* 10 ' 10"" 

frequency [HzJ 



Figure 5. Two insulator transfer functions estimated (noisy curves) and fits to 
them in the lower frequencies (thicker lines) . The inset contains data for the black 
curve — see text for details. 
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the magenta curve corresponds to a data stream about a month long and in largely 
varying thermal conditions outside the insulator; the black curve corresponds to a 
shorter time series (five days) which began with an abrupt temperature transient, 
followed by a periodic signal of a few fractions of a milli-Hertz, — see figure inset. 
These curves are relatively clean below 10""' Hz but they (noisily) tend to approach 
unity at higher frequencies, which is an indication that electronic noise is dominant in 
that band: indeed, system readouts inside and outside the insulator tend to equal each 
other (transfer function ncaring unity), while real temperature fluctuations obviously 
do not. 

The thicker lines are fits to the data in the lower frequency band: the red is the 
exact prediction of the model, while the blue is a first order filter fit, and the green 
is a second order filter fit. The last two are provided as examples that the data can 
also be adjusted by simpler models in restricted data regions yet indicating that the 
actual behaviour of the insulator follows a trend with a steeper slope towards higher 
frequencies. 

Note that the LTP MBW is above 10"'^ Hz, hence is not properly covered by 
these data. However the fact that the model predictions are followed quite well at 
low frequencies is reassuring, in that we can expect filter supression factors of 10^ and 
more in the MBW. This is simply because we of course do not expect the transfer 
function to bounce back up again at high frequencies. 

The reported facts do not constitute a full quantitative experimental test of the 
model, but they do confirm that fluctuation damping is well below the required level 
in the MBW. 

7. Conclusions 

Thermal fluctuations in the LTP must comply with very demanding requirements, 
as reflected by equation Q. Accordingly, very delicate sensors and associated 
electronics must be designed and built if meaningful temperature measurements are 
to be performed in flight. The thermal diagnostics system must then be tested on 
ground before launch. 

However, even the best laboratory conditions fall orders of magnitude short of the 
above requirement, so meaningful tests of the temperature sensing system cannot be 
tested without suitably screening the sensors from ambient temperature fluctuations. 
We have addressed how this can be accomplished by means of an insulating system 
consisting of a central metallic core surrounded by a thick layer of a very poorly 
conducting material. The latter provides good thermal insulation, while the central 
core, having a large thermal inertia, ensures stability of the sensors' environment. 
The choice of materials is flexible, so aluminium and polyurethane, which are easily 
available, have been selected. 

The appropriate sensors for our needs are temperature sensitive resistors — more 
specifically thermistors, also known as NTCs. It appears that, because these sensors 
need to be wired to external electronics, heat leakage through such wires is an effect 
which needs to be quantitatively assessed as well. We have analysed this effect, and 
concluded that thermal leakage depends strongly on the central metallic core size, and 
requires it to be somewhat large. 

Ambient temperature fiuctuations, determined by dedicated in situ measure- 
ments, are of the order of 10~^K/\/Hz at 1 mHz at our laboratory, and decrease 
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at higher frequencies within the MBW. The required stabihty conditions at the sen- 
sors, attached at the core's surface, thus need an attenuation factor of 10^, or better. 
Our analysis determines that a central aluminium core of about 13 cm of radius, 
surrounded by a concentric layer of polyurethane 15-20 cm thick, comfortably pro- 
vides the needed thermal screening which guarantees a meaningful test of the sensors' 
performance. 

It could be argued that if differential rather than absolute temperature 
measurements were performed in on-ground sensor tests, then a significant reduction 
in thermal insulation requirements would ensue. However, the OMS needs an 
absolute temperature stability of 10^^ K/\/Hz, hence an absolute temperature stability 
measurement on-ground is also needed. This is why we have chosen the insulator 
concept described in this paper. 

Our results are based on modelling. Direct experimental verification of the 
theoretical predictions in the LTP MBW is however not immediately obvious, since 
their successful implementation must result in the generation of a thermal environment 
which is more quiet than the temperature sensing system itself. Nevertheless, analysis 
of low frequency real data (i.e., below ImHz), where the insulating capabilities are 
less efficient, shows by semi-quantitative extrapolation that the insulator behaviour in 
the MBW is in practice much better than required. 

The research presented in this paper specifically applies to the LPF mission. 
But every endeavour regarding LPF has of course an ultimate motivation for LISA. 
Thermal diagnostics are no exception to the general rule and, as already pointed out, 
their design requirements do take into account that LISA will be even more exigent 
as regards thermal stability and temperature measurements. The results presented 
in this paper show that the creation of thermally very stable environments is not 
a difficult problem, and that insulator performance can be characterised with great 
accuracy by the methods presented here. In fact, the more difficult problem — and a 
major one indeed — is the design and manufacture of an electronics which be highly 
quiet down frequencies of 10"'' Hz and below, as required for LISA. We are currently 
working on these matters, and will report shortly on progress elsewhere. 
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Appendix A. Thermal insulator frequency response functions 

Here we present some mathematical details of the solution to the Fourier problem, 
equations H18|) - (|21|) . We first of all Fourier transform equations H18|) and (|21|l : 




(A.l) 



OO 




(A.2) 
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Equation (|A.1|1 can be recast in the form 

(V^ + 7^) f (x, cj) = , < r < ai (A.3) 

(V^ + 7I) f (x, Lu)=0 , ai < r < 02 (A.4) 
where r = |x|, and 

7i = -tu i— , 72 = -luj ^ A.5) 

Kl K2 

To these, matching conditions at the interface"'' and boundary conditions must 
be added: 

f{r ^ ai~0,uj) ^f{r ^ ai+0,uj) (A.6) 

df df 
Ki —{r ^ ai ~ 0, uj) = K2 —{r = ai + 0, w) (A. 7) 

or or 

f{r^a2,u})^fo{e,'f;u;) (A.8) 



Equations (|A.3|1 and (jA.4p are of the Hehiihohz kind. Their solutions are thus 
respectively given by 

Aijn{uj)ji{jir)Yi,n{0,ip) , < r < ai 

T(x,c^)=<( (A.9) 
[C(m(w) j;(72r) + Am(w) j/i(72r)] Yim(6', 1^) , ai < r < 02 

where ji and are spherical Bcssel functions |17) , 

, / 1 d V sinz , / 1 d Y cosz /. i«n 

and the coefficients A;,„(tt'), ^^(tt') and Di„i{ijj) are to be determined by 
equations HA.6|) - (|A.8p . These can be expanded as follows, respectively: 

Aim{uj) jl{liai)Yijn{0,Lp) = 

Ira 

= 53 ['^'™(^)-?''('^2ai) + Am(w)2/i(72ai)] >^im(6',<^) (A.ll) 

Im 



Im 

= K272 X! [^i™(^)ii(72ai) + A,n(^^)yK72ai)] Yi„,{0,Lp) (A. 12) 



= ^ bl„Myim{e,^) (A. 13) 



The temperature and the heat flux should be continuous across the interface. 
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Because of the completeness property of the spherical harmonics, the above 
equations completely determine the coefficients Aimi^jj), Ci,n{i^) and Di,n{u!). The 
result is 

Alr,i{uj) = (,i{uj) him{uj) , Ci™(tj) = 77; (w) bi,n{oj) , Di,n{oj) = Q{uj) bi^{uj) (A. 14) 
with 

^'(^) A ] \ ['«272 j;(72ai) yz(72ai) - K272 j/(72ai) y;(72ai)] (A. 15) 



'^'(^) [k272 j/(7iai)y/'(72ai) - ki7i J; (7iai) yi(72ai)] (A.16) 

^'(^) = A ] \ [«i7i ji(72ai) j;(7iai) - K272 j;'(72ai) j;(7iai)] (A.17) 



and 



Ai{uj) = Ki7i j,(7iai) [j; (7201) y; (7202) -j! (7202) 2/; (7201)] + 

+ K272i/(7iai) [i/(72a2)2/i(72ai) - j[(72ai)y;(72a2)] (A.18) 

When the above results, equations IjA.lSp through (|A.18p . are inserted back into 
equation (jA.Op the result stated in equation H22|l in the main text obtains, i.e., 

f(x,w) =^ i?i„(x,a;)6i„(a;) (A.19) 

Im 

where 

„ , , r aHjK7iO>1m(^,^) , 0<r<ai 
i7i,„(x, w) = < (A.20) 
t Ji(72'') +0(^^)2/i(72?')] yiml^*, <^) , ax<r <a2 

For monopole only boundary conditions, equation l|28(l . the transfer function is 

rr. N [ ^o(^^)jo(7iO , 0<r<ai 
H{r,uj) = < (A. 21) 

t m{^)jo{l2r) + Co(w) yo{l2r) , ai<r <a2 
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